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A flavour symmetry based on Z4 is analysed in the context of SU(5) Grand Unification with the 
standard fermionic content plus three right-handed neutrinos. The role of Z4 is to forbid some 
Yukawa couplings of up- and down-quarks to Higgs scalars such that the quark mass matrices 
M u , Md have Nearest-Neighbour-Interaction (NNI) structure, once they are generated through the 
electroweak symmetry breaking. It turns out in this framework that Z4 is indeed the minimal 
discrete symmetry and its implementation requires the introduction of at least two Higgs quintets, 
which leads to a two Higgs doublet model at low energy scale. Due to the SU(5) unification, it is 
shown that the charged lepton mass matrix develops also NNI form. However, the effective neutrino 
mass matrix exhibits a non parallel pattern, in the framework of the type-I seesaw mechanism. 
Analysing all possible zero textures allowed by gauge- horizontal symmetry SU(5) x Z4, it is seen 
that only two patterns are in agreement with the leptonic experimental data and they could be 
further distinguished by the light neutrino mass spectrum hierarchy. It is also demonstrated that 
Z4 freezes out the possibility of proton decay through exchange of colour Higgs triplets at tree-level. 
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I. INTRODUCTION 

Grand Unified Theories (GUTs) are beautiful attempts 
beyond the Standard Model (SM) for understanding the 
observed quark and lepton masses and their mixings, the 
so-called "flavour puzzle" . This is indeed corroborated by 
the fact that the running gauge couplings when evolved 
to very large energy scales, typically 10 14 ^ 16 GeV, seem 
to unify to a unique coupling. The simplest GUT model 
which accommodates the SM gauge and fermion fields in 
a few multiplets is based on thegroup SU(5), proposed in 
1974 by Georgi and Glashow [lj. Some small extensions 
of the Georgi-Glashow model are still viable today 
even in its supersymmetric version 

Generically in GUT models, not only the SM gauge 
couplings do unify but also the SM fermions are unified 
in a small number of large multiplets that lead to new 
phenomenological signatures. An important signature of 
most GUTs is the prediction for proton decay , which 
has not yet been observed and then severely constrains 
the GUT models. The fact that the quarks and leptons 
are tied together in GUT multiplets is not enough to fully 
determine the properties of their observed masses and 
mixings. However, the GUT relations among quark and 
lepton Yukawa matrices are an excellent starting point 
for building a flavour symmetry. Thus, if one requires a 
flavour symmetry to enforce a particular pattern in the 
up- and down-quark Yukawa couplings, this engenders 
physical consequences in the leptonic sector. 



During the last decades, a huge number of flavour 
symmetries with different purposes have been extensively 
presented in the literature (in the context of SM, Grand 
Unification, etc.). The simplest and attractive possibility 
is to assume the vanishing of some Yukawa matrix ele- 
ments ( "texture zeroes" ) by the requirement of a discrete 
symmetry [12H17I ]. such that it would naturally lead to 
the flavour mixing angles be expressed in terms of mass 
ratios. The converse is not necessarily true, since one 
can obtain zeroes in the Yukawa matrices fl8H20| just 
by performing some set of transformations (weak basis 
transformations) leaving the gauge sector diagonal. For 
instance, it is remarkable in the SM that one can al- 
ways go to a weak basis where both up- and down-quark 
mass matrices M u , Md have simultaneously the form or 
a "parallel structure" , 
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known as the Nearest-Neighbour-Interaction (NNI) ba- 
sis [2l[ . Being the matrix form in Eq. ([!} for both quark 
sectors just a weak basis, no zero in the NNI matrices has 
physical meaning. The NNI basis is closely connected to 
the Fritzsch ansatz [22I - I24I ] . which further assumes the 
Hermiticity condition on the NNI quark mass matrices 
M„, Md- Through a simple choice of weak basis trans- 
formation it is always possible to make M u , Md Hermi- 
tian, but their structures are no longer of the NNI form. 
Assuming that M u , Md are in the NNI basis, it has been 
shown in Ref. 25] that the experimental data are still 
in agreement with relatively small deviations from Her- 
miticity, at the 20% level. This procedure was also ex- 
tended to the leptonic sector in Ref. [2t| . 
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Furthermore, it was also shown in Ref. [HJ that it is 
possible to attain the up- and down-quark mass matrices 
M u , Md with NNI structure through the implementation 
of an Abelian discrete flavour symmetry in the context of 
the two Higgs doublet model (2HDM). In that context, 
the minimal realisation is the group Z4. In a general 
2HDM, a NNI form for each Yukawa coupling matrices 
cannot be a weak basis choice. Indeed, the requirement 
of the Z4-symmetry does imply restrictions on the scalar 
couplings to the quarks, although one gets no impact on 
the quark masses and the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix [27|, H|. 

The purpose of this article is to study whether it is pos- 
sible to construct a Z4 flavour symmetry, similar to the 
one implemented in Ref. [25|, that leads to quark mass 
matrices M u , Md in the NNI form in the context of SU(5) 
Grand Unification with the usual fermionic content. In 
addition, since SU(5) implies relations among quarks and 
leptons, we also explore the physical consequences of the 
Z4 flavour symmetry on the leptonic sector for the case 
where three ri ght- handed neutrinos are added, being the 
type-I seesaw [29|-|32j the mechanism responsible for the 
light neutrinos to acquire Majorana masses. 

This article is organised as follows. In Section [TT| we 
introduce the SU(5) x Z4 model. Next, in Section [TTT1 we 
analyse the different channels of proton decay as well as 
the issues of unification in this model. Then in Section lTVl 
we discuss in detail the form of the leptonic mass matrix 
textures provided by the flavour symmetry. In Section IVl 
we present our numerical analysis on the leptonic sector 
as the result of the restrictions due to Z4 . The zero tex- 
tures obtained for the effective neutrino mass matrix are 
then confronted with leptonic observable data. Finally, 
our conclusions are drawn in Section fVTl 



II. THE MODEL 

Following Ref. (25|, we build an Abelian discrete 
flavour symmetry within a SU(5) GUT model which 
yields at low energies to a NNI form for both up- and 
down-quark mass matrices. We choose the flavour sym- 
metry to be Abelian, because it is the simplest way to for- 
bid some Yukawa couplings, so that texture zeroes appear 
naturally in the mass matrices. The flavour symmetry is 
also chosen to be discrete in order to avoid the presence 
of Nambu-Goldstone bosons. To simplify our search for 
a minimal flavour symmetry realisation on the full La- 
grangian, we consider only the case where the flavour 
group belongs to the Z„ family. Thus, each fermionic or 
Higgs multiplet, R, transforms as 



R 



(2) 



where the charges Q(R) € Z n . 

The particle content of our GUT model is a small ex- 
tension of the original SU(5) model proposed by Georgi 
and Glashow [l[ in 1974. It contains three generations of 



10, 5* fermionic multiplets, which accommodate the left- 
handed fcrmions of the SM, Q it it?, d\, L i7 e?, as follows 



10* = (Q l ,u c l ,ef) , 5* = (Li,dl) 



(3) 



where i = 1,2,3 stands for the generation index. Fur- 
thermore, we introduce three right-handed neutrinos is? 
(in the left-handed picture), singlets under SU(5), as the 
simplest way to generate the light neutrino masses needed 
to explain the observed neutrino oscillation data. Being 
not constrained by any gauge symmetry, the singlets vf 
can have a Majorana mass term and a Dirac Yukawa term 
mixing with the leptonic doublets. After the electroweak 
symmetry breaking, if one takes the Majorana mass term 
close to the GUT scale, so much larger than the Dirac 
mass term, light neutrino masses can be generated - the 
type-I seesaw mechanism [29-32]. 

The Higgs sector of the model consists of an adjoint 
Higgs multiplet, £(24), chargeless under Z„, and two 
quintets Hi(5), i?2(5), with different Z„ charges </>i, <^2, 
respectively. Under these charge assignments the full 
scalar potential V reads as 



U = -i M 2 Tr(E 2 ) + iaTr(i: 



3\ 



h? [Tr(£ 2 )] 2 + ^Tr(£ 4 ) 

H\ Qmi + aiE + An Tr(£ 2 ) + A 12 £ 2 ) Hi 



H 2 ( o A*2 + «2E + A21 Tr(£ 2 ) + A 22 £ 2 ) H 2 



1 

2' 

A 1 |-ff 1 | 4 + A 2 | J ff 2 | 4 + A 3 |ffi| 2 l-ff2| 2 



X 4 [H{H 2 H^H 1 



where the parameters n, a and b have mass dimensions, 
while A is dimensionless. Note that the self-potential 
terms for the adjoint field £ is as general as in the mini- 
mal SU(5) model. 

The adjoint field £ breaks spontaneously the SU(5) 
gauge group to the SM group (SU(3) C x SU(2) L x U(1) Y ), 
through the vacuum expectation value (VEV), 



<£) =adiag(2,2,2,-3,-3), 



(5) 



provided that a is the following solution that minimises 
the scalar potential [33T - [36| given in Eq. (Q}, 



1 + y/l + 4^ (6077 + 7) 



2A 



6O77 + 7 



(6) 



where 77 = b 2 /X and £ = A/1 2 /a 2 requiring that 77 > 
—7/60 and A > 0. The parameters 77 and £ are enough to 
allow whether the VEV from Eq. ([5]) corresponds to an 
absolute minimum of the potential. Thus, the parameters 
of the potential have to be properly chosen to guarantee 
the SM group in the broken phase [35[ with the natural 
value for a lying around the unification scale A. The 
VEV given in Eq. also splits the adjoint Higgs field £ 
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in its components E3 (weak isospin triplet), Ss (colour 
octet) and S 2 4 (singlet), which become massive. 

The Higgs quintets Hi, H 2 are introduced to break 
the SM gauge group down to SU(3) C x U(l) c . m . and also 
generate the fermion masses via the Yukawa interactions 
at the electroweak scale. The quintets Hi, H2 are split 
by the VEV (£) into the Higgs doublets $1, $ 2 and the 
Higgs colour-triplets Ti, T 2 , respectively. Thus, the SM 
group is then broken through the VEVs vi , v 2 of the re- 
spective Higgs doublets $i,$ 2 , verifying 



+ M 2 = (V2Gf) 1 = (246.2 GeV) 2 , (7) 



where Gp is the Fermi constant. 

Furthermore, one has to avoid rapid proton decay me- 
diated by the Higgs colour-triplets Ti, T 2 , which can be 
solved by fine-tuning the parameters of the Higgs poten- 
tial, O (f ) ~ io-( 12 - 13 ) - the so-called doublet-triplet 
splitting problem. This fine-tuning can be re-expressed 
as new constraints on the mass parameters n\ and fx 2 , 



6<7 (ai — lOerAn — 3ctAi 2 ) 
6a (a 2 - 10crA 2 i - 3ctA 22 ) 



(8a) 
(8b) 



such that the Higgs colour triplets Ti, T 2 have masses: 



= 5a(ai - crAi 2 ) , 
m% = 5a(a 2 - crA 22 ) . 



(9a) 
(9b) 



Thus, once the heavy Higgs colour-triplets Ti, T 2 are in- 
tegrated out, the model is just a two Higgs doublet model 
(2HDM). Since the adjoint Higgs multiplet carries no Z„ 
charge, the obtained 2HDM automatically preserves the 
flavour symmetry in higher orders of perturbation the- 
ory, provided that no Nambu-Goldstone boson ap pea rs 
at tree- level due to an accidental global symmetry [33]. 

In order to fully determine the Z„ charges for the 
fermions and the Higgs quintets such that the quark mass 
matrices M u , Md have NNI form, one needs to analyse 
the Yukawa interactions [25]. The most general Yukawa 
Lagrangian reads as 

= \ (ri) tf 10, 10, Hi + \ (r 2 ^. 10, 10, H 2 



V2 (rj) . . 10, 5; hi + V2 (t 2 ).. io, 5* h; 
(r^).. 5*^Hi + (r 2 D )„ 5*^h 2 



(10) 



ij « 3 



where r*' and Mr are symmetric complex matrices, 



while T d ,2 D are just general complex matrices. The quark 
mass matrices M u , Md are then given by 

M u =vi Ti + v 2 T 2 U , (11a) 
M d =v*iT l d + v* 2 T 2 , (lib) 

and their zeroes are directly settled by the Yukawa matri- 
ces r^' 2 . It is clear from Eq. (|lla|) that the mass matrix 
M u is symmetric. 



In order to guarantee that the (33)-element of up-quark 
matrix M u does not vanish we set the charge 2 as 



-2<? 3 



(12) 



where q 3 = Q(103). This choice automatically fixes the 
Z„ charges of the multiplets 5*, 10^ as a function of 0i 
and (73 as 



Q(10j) = (3q 3 + <j)i, -q 3 - </>!, q 3 ) , 
Q(5*) = (<? 3 + 20i, -3<? 3 , -q 3 + 0i) . 



(13) 



The Z„ charges of the right-handed neutrinos, vf (singlets 
under SU(5)) have no restrictions and are taken as free 
parameters, 



(14) 



It is in fact not enough to have the charge assignments 
given in Eq. (TT3|) in order to guarantee the NNI structure 
for the mass matrices M u , Md, one has in addition to pre- 
serve the NNI-zero entries, which implies that one should 
forbid the quark bilinears corresponding to the NNI-zero 
entries to couple to Higgs doublets. The analysis of the 
quark bilinears can be made directly from the 10^10^ and 
10^5* bilinears. By taking into account Eq. (J3|) one can 
derive the 10; 10, bilinear charges as 



'6g 3 + 20i 2q 3 Aq 3 + cpi 
2q 3 -2^1 - 2g 3 -0i 
4g 3 + 0i -0i 2g 3 



(15) 



and the 10i5* bilinear charges as 



'4^ + 301 0i 20i + 2^ 

0i -4g 3 - 0i -2^3 
20i + 2q 3 -2q 3 0i 



(16) 



One sees immediately from Eq. (fTS"]) , that in the case of 
0i = 2 one gets that both Higgs doublets $1, $ 2 can 
couple to the bilinear IO1IO1, thus destroying the NNI 
form for M u . According to the Eqs. (|15[) and (HHJ), the 
lowest order group of Z n -type which res pec ts the zero- 
entries in the NNI form is Z4, as in Ref. [25|]. The sym- 
metry group Z4 is indeed the minimal discrete flavour 
symmetry that gives rise to NNI structure for M u , Md, 
since the other candidate, Z 2 x Z 2 , having all its non- 
trivial elements of order 2, is not viable. 

The zero-entries in the Yukawa coupling matrices F^ 
are unequivocally determined by Eq. (|13[) . provided that 
the bilinear entries given in Eqs. (|T5|) and (jTHJ) are taken 
correctly into account, as follows: 



'0 

ri = I b u 

,0 b u 



r 1 - 

1 d — 



a d 
a'd 
c dl 




(17a) 



(17b) 
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which when applying Eqs. (fTTj) leads to mass matrices 
M u ,M d with a NNI form. We point out that the NNI 
structures thus obtained have a different nature than in 
the pure SM case. As we have emphasised in the intro- 
duction, in the SM the NNI structure for M u , Md is just a 
choice of weak basis. Instead in our model, the NNI form 
for M u , Md arises as the requirement of a Z4 symmetry, 
i.e. the form of the Yukawa couplings P*'^ are dictated 
by Z4 in Eqs. (fTTl) . and therefore it is not a weak basis 
choice. In fact, the Z4 do imply new restrictions on the 
scalar couplings to quarks. 

In order to complete the construction of the flavour 
symmetry Z4, one has to ensure in addition that the 
scalar potential does not acquire an extra accidental 
global symmetry which, upon spontaneous electroweak 
symmetry breaking, would lead to a massless Nambu- 
Goldstone boson at tree level (38[- It is indeed true 
that the most general scalar potential invariant under 
SU(5) x Z 4 , given in Eq. possesses an extra U(l) 
global invariant transformation on the quintet fields. 
Analogously of what was described in (25[, this problem 
can be cured by soft-breaking the Z4 symmetry through 
the introduction of a term like 



(18) 



in the scalar potential in Eq. Notice that the soft- 
breaking term /112 is not sufficient to break CP spon- 
taneously in this model, so that CP violation arises 
simply from complex Yukawa couplings, leading to the 
Kobayashi-Maskawa mechanism, see KM in Ref. [28[ . 

An elegant way to avoid the term given in Eq. (|18[) and 
at the same time prevent the existence of the global U(l) 
symmetry is simply by adding a complex SU(5) singlet 
Higgs field S which transforms non-trivially under Z4. 
The most general potential involving the singlet field S 
reads as 



V S 



H\ (M12 + A' 12 £)ff 2 S + H.c 



1 



+ X S \S\ A + X' S (S 4 + H.c). 



(19) 



Note that the last term in the Eq. (|19j) . allowed by Z4, ex- 
plicitly prevents the appearance of Goldstone bosons that 
could result from an accidental global U(l). In the case 
of the VEV of the complex singlet field S is of the order 
of the GUT scale, it is easy to verify that arg((5)) = 
or 7r/4, but it does not break the CP symmetry sponta- 
neously. 

We are now able to derive the texture zeroes in the 
leptonic sector from the constructed SU(5) x Z4 model. 
Since in this model the relation 



M ' = Mj 



(20) 



holds at the GUT scale and the zeroes of the NNI struc- 
ture are placed symmetrically, one easily sees that the 
flavour symmetry Z4 also leads to a charged lepton mass 



matrix M e with NNI form like for quark mass matrices 
M u , Md- On the other hand, for the neutrino sector we 
do not expect a parallel structure neither for the Dirac 
mass matrix, mo, 



vi T 1 D + v 2 T 



D ' 



(21) 



nor for the Majorana mass matrix Mr, since the Z4 
charges of the right-handed neutrino are taken arbitrary. 

It is well-known that the GUT relation given by 
Eq. (|20[) is not in fact compatible with the down-type 
quark and the charged lepton mass hierarchies observed 
at low energies (badly violated for the first and second 
generations) [39j |. In the following subsections, we shall 
roughly sketch two small extensions of the model that 
have the aim to modify properly the GUT relation given 
in Eq. (|20[) without changing the zeroes in the quark 
mass matrices M u , Md and the charged lepton mass ma- 
trix M e . 



A. "Consistent" SU(5) 

A possible way to modify correctly the relation given 
by Eq. (f2"U]) without adding new representations is by 
considering some non-renormalisable higher dimensional 
operators [7|,l9|, 0(1/A'), due to physics above the GUT 
scale, A' ^> A. The natural scale for A' is to be two or 
three orders of magnitude greater than the GUT scale. 
Higher dimensional operators involving the adjoint field 
£ of the type 



n=l,2 



V2 
A' 



(An)* 



H* na 10f ££5* c) 



(22) 



with the indices a, 6, c — 1, ■ ■ • ,5 and i,j = 1,2, 3, would 
contribute to the relation in Eq. (|2TJ|) as 

M d - Mj = Ai + v* 2 A 2 ) . (23) 

The complex matrices Ai and A 2 can account for the dis- 
crepancies between M d and M e . The quark and charged 
lepton mass matrices M u ,M d ,M e remain in the NNI 

form, since the adjoint field £ is trivial under Z4 and 

1 2 

contributes to the quark Yukawa matrices r u ' d through 
dimension-five operators when it acquires VEV. The only 
higher dimensional operators that could spoil the NNI 
structure on the quark mass matrices are of dimension- 
six, e.g. 



A 



A 



-2 10 2 10 2 ^ H* H 2 , 



(24) 



and therefore very much suppressed, 0(v 2 /A' 2 ) . More- 
over, the presence of dimension-five operators in the 
Higgs potential given in Eq. (j4} can contribute to the 
splitting between the masses of the Higgs multiplets £3 
and £g by several orders of magnitude. 
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B. Adjoint SU(5) 

The other alternative consists in maintaining the full 
Lagrangian renormalisable just by requiring a 45 dimen- 
sional Higgs scalar 6], 7^(45), instead of the quintet 
#2- The field representation "H"* 3 satisfies the relations: 

H"? = -H^ a and ^Li U a J = . Thus, the potential 
given in Eq. is now modified by the following terms: 



f 



^ + Aai Tr(S^) + a 2 S + A, 2 ^ J j H 

+ A 2 (H t H) 2 + A 3 |i?i| 2 (H t H) 



A 4 [Hlm^Hx 



that substitute the terms involving the quintet H2. For 
the sake of simplicity of the notation, the SU(5) invariant 
contractions involving % were not explicitly written. The 
mismatch between Md and M e is explained as 



M d - Ml = 



d u i5 > 



(26) 



where W45 is the strength of the vacuum expectation value 
of the field T-L, assuming 



/35\ 



V45 



(27) 



It is clear that if the Higgs multiplet H has the same Z4 
charge as the one assigned for H2 in Eq. (|T2"j) , one recovers 
the NNI form for the quark mass matrices M u , Md and 
the charged lepton mass matrix M e . 



III. PROTON DECAY AND UNIFICATION 

In this section we analyse the proton decay in the con- 
structed SU(5) x Z4 model. In this model, the proton 
decays through the exchange of the heavy lepto-quark 
gauge bosons X, Y or the colour Higgs triplets T\,T2- 
The experimental limits on the proton decay rate severely 
constrain the masses of such heavy states that we shall 
assume of the order of the unification scale A, since in 
this scenario the proton decay rate is inversely propor- 
tional to the mass square of the heavy states. On the 
other hand, the unification scale is by definition the scale 
where the running gauge couplings measured at the scale 
M z = 91.1876±0.0021 GeV [H do unify. Thus, the lim- 
its on the proton decay rate have to be confronted with 
the parameters that govern the evolution of the gauge 
couplings. 

The twelve lepto-quark gauge bosons X, Y (compo- 
nents of a colour weak isospin doublet) arise from the 
adjoint 24 representation that also contains the twelve 
gauge bosons of the SM. The gauge bosons X, Y become 
massive through the Higgs mechanism with a common 
mass, My, 



My 



25 



2 J2. 



-9u° 



(28) 



where gu is the unified gauge coupling. To suppress the 
X, Y boson proton decay channels, one has necessarily 
that My ^> m p (the proton mass) which then leads to 
an approximate four fermion interaction (dimension-six 
operators) proportional to 1 /M v and the unified coupling 
ay = /Ait . In this approximation, the proton decay 
width can be estimated as 1411: 



2 m P 

M^ 



(29) 



Making use of the most restrictive constraints on the par- 
(25) tial proton lifetime r(p 7r°e+) > 8.2 x 10 33 years jioj . 
one can derive a rough lower bound for the X, Y mass 
scale My, 



M v > (4.0-5.1) x 10 i5 GeV, 



(30) 



which corresponds a range of the unified gauge coupling 
a" 1 fa 25 — 40, as suggested by performing the renor- 
malisation group evolution of the gauge couplings (see 
details in Appendix \K§ . Since we assume for the unifi- 
cation scale A ~ My, the constraint given by Eq. ([50)1 
determines the scale where t he g auge couplings should 
unify (for a recent review see 11|). 

Usually in non-supersymmetric scenarios, the proton 
decay through the exchange of Higgs colour triplets T\ , T2 
is very suppressed. Being these decay modes also de- 
scribed by dimension-six operators, their suppression is 
proportional to products of Yukawa couplings, which are 
much smaller than the gauge couplings in the X, Y boson 
exchange. In fact, the contribution of these dimension-six 
operators vanishes at tree-level when the Z4 symmetry is 
exact. The dimension-six operators contributing to the 
proton decay via the colour triplet exchange are given at 
tree-level by: 



E 

,1=1.2 



kl 



(QiQjXQkLi) + (u?ep«d?) 



(31) 

It is then clear from the pattern of the Yukawa coupling 
matrices and T d given in Eqs. (fl~7|) that the only pos- 
sible non- vanishing contribution of the dimension-six op- 
erators given in Eq. (|31[) involve necessarily fermions of 
the third generation. One concludes that at tree-level the 
proton does not decay through the four-fermion interac- 
tions described by the operators given in Eq. (|31[) . 

If one soft-breaks the Z4 symmetry of the potential 
with a /ii2-term as the one written in Eq. (|18[) , such 
term mixes already at tree-level the heavy states T\, T2 
and therefore induces proton decay through dimension- 
six operators proportional to (r 2 )^. (r^)^, thus involv- 
ing fermions of the first and the second generations. This 
can be avoided when a singlet scalar S, charged under Z4, 
is introduced in the potential as shown in Eq. (|19l) and in 
this case proton decay via Ti , T 2 exchange vanishes once 
more at tree-level. 

In what concerns unification, it is widely established 
that the running gauge couplings do not unify in the 
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FIG. 1. The plot of the gauge coupling evolution, ot^\ 3 , at 
two-loop level for the choice of Ms 3 = 500 GeV. The uni- 
fication then occurs for A = 1.3 X 10 14 GeV and Me. = 
3.2 x 10 11 GeV. The dashed line corresponds to the running 
of the gauge couplings without considering the intermediate 
scales Me, and Me s • 



context of the SM even in the presence of an extra Higgs 
doublet [ll|. However, if one considers the splitting be- 
tween the masses of the multiplets £3 and Eg, it turns 
out possible to unify the gauge couplings at two-loop level 
even without taking into account threshold effects. We 
looked for unification assuming the masses of the bosons 
X, Y, T\ and T2 at the unification scale A and allowing a 
splitting between the masses M% 3 , M^ s ■ We also assume 
the masses of the Higgs doublets $1, $2 lying around the 
electroweak scale. 

From our numerics, we have found a small variation 
for the unification scale A, 

1.3 x 10 14 GeV < A < 2.4 x 10 14 GeV , (32) 



and for the masses Ms 3 and Ms 8 , 

M z < Ms 3 < 1.8 x 10 4 GeV, 
5.4 x 10 11 GeV < M E . < 1.3 x 10 14 GeV , 



(33a) 
(33b) 



which corresponds to a mass difference of the order 
Ms 8 /Me 3 = 0(1O 9 - 10 ). These results do not get im- 
proved even when a large splitting between the two Higgs 
doublets is considered. The details concerning the equa- 
tions of gauge coupling evolution are sketched in Ap- 
pendix [5] For illustration, in Fig. [T] we plot the gauge 
coupling evolution at two-loop level for Mj 3 = 500 GeV, 
which fixes A = 1.9 x 10 14 GeV and M Sg = 3.2 x 10 12 
GeV in order to achieve unification. 

The unification of the gauge couplings suffers from two 
potential problems. First, the scale A in Eq. (|32|) seems to 
be lower than what is required by the lower bound given 
in Eq. (|30p. Second, the mass splitting between Mj 3 and 
Ms s is unnaturally large. This may suggest the necessity 



of introducing extra multiplets in order to relax such con- 
straints, e.g. see Rcfs. [2, 6]. It is interesting to remark 
that the alternative presented in subsection III Bl can suc- 
cessfully adjust the running of SM gauge couplings, since 
the broken components of the 45 Higgs representations 
have the correct quantum numbers [6| for the unification 
of the gauge couplings. 



IV. QUARK AND LEPTON MASS MATRICES 

In Section [Til we have derived a Z4 flavour symmetry 
that induces the quark mass matrices M u , M4 to have a 
NNI structure in the framework of SU(5) Grand Unifica- 
tion with minimal fermion content. Due to the fact that 
quarks and leptons are tied together in fermionic multi- 
plets, the flavour symmetry imposed at the Lagrangian 
also restricts the form of the charged lepton mass matrix 
M e , which develops a NNI structure. Moreover, the ex- 
tensions considered in subsections III Al and III Bl not only 
lead to Md ^ Mj but also M u is no longer symmet- 
ric @, H, El|. Hence, in what follows we shall assume 
arbitrary NNI mass matrices M u , Md and M e , where the 
asymmetry among the elements A u _d,e and A' u d e , B Ui d, e 
and B' u d can be measured by the deviation parameters 

u,d,e 

e„ V as 



A' 




A u ,d,e 


A 1 


+ 




B' 

u,d,e 




B A 

u,d,e 


T3I 

u,d,e 


+ 


B A 



(34a) 



(34b) 



where the parameters A e , A' e , B e and B' e follow the con- 
vention made in Eq. ([T]) for the quark sector. Global 
measurements of the asymmetry in the quark, e q , and 
charged lepton, e e , sectors are defined as 



(et) 2 



(35a) 
(35b) 



It is interesting to see that in the limit when e q = e e = 
one recovers the Fritzsch ansatz [22Tj24] [ . 

Concerning the neutrino sector, since neutrino charges 
are taken as free parameters in Eq. (IT4|) . one has to clas- 
sify all viable textures for the mass matrices m u and Mu 
by scanning all combinations of the Z4 neutrino charges. 
Hence, no NNI form is expected for the mass matrices 
mo and Mr and their resulting textures have to be con- 
fronted with the neutrino experimental data. We start 
our scanning by deriving the allowed range for the Higgs 
doublet charges <f>\, 4>2- As dictated by Eq. (|12p. the 
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charge of the Higgs doublet $2 can only take two val- 
ues: (j>2 — or 2 . On the other hand, the charge of the 
Higgs doublet $1 has to be odd (0i — 1 or 3 ), which can 
be clearly seen, for instance, by noting that the entries 
(2,2) and (3,3) of the bilinear given in Eq. (|15|) would 
be equal when <pi is even. Once the Z4 charges of the 
right-handed neutrino fields, i>i, are fixed, one can im- 
mediately determine the pattern of the effective neutrino 
mass matrix m v . From the structure of the Lagrangian 
in Eq. (I10p one can derive the effective neutrino mass 
matrix m„, which is given by the usual standard type-I 
seesaw formula |23 - l32l |: 

m v — —mo Mr 1 m J) ! (36) 

to an excellent approximation. Since the symmetric Ma- 
jorana mass matrix Mr is directly introduced at the La- 
grangian level, its pattern is determined by the charges 
V{ . Thus the Z4 charges of bilincars uf Uj are given by 

2v 1 v 1 + v 2 v\ + v 3 \ 
v 2 + v\ 2v 2 v 2 + v 3 \ . (37) 
v 3 + V\ v 3 + v 2 2v 3 J 

The texture zeroes in the Dirac neutrino mass matrix m d 
are then obtained thereby computing the charges of the 
bilinears 5* , 

q 3 + 20i + vi q 3 + 20i + v 2 q 3 + 20i + 2/3 \ 
-3q 3 + z^i -3q 3 + v 2 -3(?3 + v 3 (38) 

-?3 + 01+^1 -<?3 + 01 + ^2 ~qz + 01+^3/ 

and verifying their couplings to the Higgs doublets. Fi- 
nally, one is then able to compute the texture zeroes in 
the effective neutrino mass matrix, m„, just by applying 
the seesaw formula given in Eq. (|3"6]) . Thus, by spanning 
all allowed values for the charges 0i, q 3 and Vi, one can 
draw all possible zero textures allowed by the symmetry 
Z 4 . We sketched in Table Q] all possible effective neutrino 
mass matrices m„ obtained as a function of the param- 
eters 0i, q 3 and The notation used to represent all 
textures in Table U reflects the fact that each of them is 
related to one of the following three classes of textures: 

/O A 0\ [OA 0\ (A 0\ 

I) I A B I , II) [A B C , III) B C , 
\0 B C) \0 C D) \0 C Dj 

(39) 

through a permutation matrix P g of the form 

m„ <9) = P g m v Pj , g G S 3 , (40) 

where {P g } are the six 3x3 permutation matrices iso- 
morphic to the symmetric group S 3 . This identification 
is also useful in simplifying the diagonalisation proce- 
dure. It is clear from the seesaw formula given in Eq. ([36]) 
that any permutation among the right-handed neutrino 
Z4 charges does not change the pattern of the effective 
neutrino mass matrix m v . This is indeed the reason why 



TABLE I. The obtained textures zeroes for the effective neu- 
trino mass matrices in the context of SU(5) x Z4 symmetry 
with two-Higgs doublets. 







v= (0,1,3) 


v= (1,2,3) 


Vi € {0, 2} 







1(132) 


11(12) 


HI ( i 2) 


01 = i 


1 


1(13) 


II 


III 




2 


11(12) 


1(132) 


111(12) 




3 


II 


1(13) 


III 







1(132) 


11(12) 


111(12) 


01 =3 


1 


II 


1(13) 


III 




2 


11(12) 


1(132) 


111(12) 




3 


1(13) 


II 


III 



we have denoted the charges Ui in Table U just by one 
ordered 3-tuple. 

Unlike what happens on the quark mass matrix pair 
M u , Aid, it is remarkable to verify that the seesaw for- 
mula in Eq. ()36[) together with the allowed 0i, q 3 and vi 
charges lead necessarily to a non-parallel structure in the 
leptonic mass matrix pair M e , m„ . No NNI form was 
found for the effective neutrino mass matrix m„. Never- 
theless, in the case where neutrinos are Dirac fermions, 
a parallel structure with NNI form in the leptonic sec- 
tor is indeed possible, if one also requires the flavour 
symmetry to forbid the appearance of the Majorana 
right-handed mass matrix, Mr. The minimal discrete 
group realisation is Z7 with, for example, the following 
charge assignments: (0i,02) = (1,0), Q(10) = (2,4,3), 
Q(5*) = (3,5,4) and v = (1,3,2). In this example, the 
NNI mass matrix mo encodes all neutrino masses and 
mixings and it was shown in Ref. (26| that the mass ma- 
trix pair M e , m,£, with parallel structure can accommo- 
date the leptonic experimental data. 

Analysing carefully the zero textures obtained in Ta- 
ble U some comments are in order. If at least two 
right-handed neutrinos have the same Z4 charge one re- 
alises that one generation decouples from the others (i.e., 
Texture-III and its permutation (12)), it corresponds 
to the right-handed neutrino charges having only values 
vi = or 2. Due to the fact that the charged lepton mass 
matrix is in the NNI form, it is rather easy to see that 
having a generation which decouples in the neutrino sec- 
tor is not phenomenologically viable, since it is impossible 
to obtain large mixing angles. In a similar way Class-I, 
where the effective neutrino mass matrix is a permutation 
of the NNI matrix, is not viable too since it leads to small 
mixing angles. Also by counting the number of indepen- 
dent free parameters of the mass matrix set M e , rriy in 
Class-I one gets a total of ten parameters which have to 
account for the twelve low energy lepton observables (six 
lepton masses, three mixing angles and three phases). 
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One is then left with only two zero textures of Class-II, 
namely Texture-II and -II( 12 ) . The confrontation of such 
textures with the neutrino data is explored in the next 
section. 



V. NUMERICAL ANALYSIS OF NEUTRINO 
MASS MATRIX WITHIN THE CLASS-II 

In this section we analyse the phenomenological conse- 
quences of each effective neutrino mass matrix belonging 
to the Class-II. As discussed in the previous section only 
two zero textures II and H(i2) from Table HI need to be 
confronted with the observable neutrino data. Without 
loss of generality one can write the charged lepton mass 
matrix, M e , and the effective neutrino mass matrices, 



in 



(</) 



of Class-II as: 
/ 



M K = k! 











MI-ZD 







C R 



Pn 



( A v 

Ay By Cy 

\ C v D v 



P 



(41a) 



(41b) 



where g — e or (12) according to Table HI the con- 
stants A e , B e ,A Vl B„, C e ,i/, D v are taken real and posi- 
tive. The diagonal phase matrix K e can be parameterised 
as 



diag(e lKl ,e lK2 ,l). 



(42) 



and the phase tp in Eq. (|41bj) cannot be absorbed by 
any field redefinition. In the case of Texture-II( 12 ) the 
permutation matrix P(i2) is simply given by 





(o 


1 








1 







(43) 




U- 







I 



Although the number of the parameters encoded in the 
pair M e , m v is twelve as the number of independent phys- 
ical parameters experimentally observed at low energy, 
the zero pattern exhibited in Eqs. (|4ip does imply new 
constraints among the independent physical parameters, 
as it will be shown. 

In order to extract the mixings angles and the CP 
phases encoded in the charged lepton and the light neu- 
trino mass matrices in Eqs. (|4ip . one needs to evalu- 
ate the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) ma- 
trix [42h44| , U. It is then useful to introduce the Hermi- 
tian mass matrix H e defined by: 

H e = M e Ml . (44) 
From Eq. (|41a|) one deduces that H e is a real Hermitian 
mass matrix, which can be diagonalised by a real and 
orthogonal matrix, O e , in the following way: 



Oj H e O e = diag(mjj, m 2 m\) 



(45) 



To diagonalise the effective neutrino mass matrix m£? it 
is practical to define the complex matrix m[J as 



where m[J is diagonalised by U v as, 

Uj ml U v = diag(mi, m 2 , m 3 ) , 

and m% are the positive light neutrino masses, 
the PMNS matrix U is given by 

U = Oj K\ P q U v . 



(46) 

(47) 
Finally, 

(48) 



It is useful to re-express the unitary matrix U in terms 
of the standard parameterisation, which has the property 
to better express the neutrino observed data in a more 
clear and uniform way, 



U 



C12C13 

-S12C23 

S12S23 - Ci 2 C23Sl 3 e 



Cl2S23Sl3e 

iS 



C12C23 " 
-C12S23 



S12C13 
Sl2S23Sl3e l<5 
- Sl2C 2 3Sl3e l 




(49) 



where Si 



6 is a Dirac CP 



violation phase and a,\, ai are Majorana CP viola- 
tion phases. From the neutrino oscillation experiments 
one infers [45j the neutrino mass squared differences 



A? 



l 21 ) 



Ar 



'31 



(Am 



m~) as well as the mixing 



angles 812 , 823 and #13 that are shown in Table HB 

In addition to oscillation parameters, we have also con- 
sidered in our numerical program three other constraints: 



the effective Majorana mass, m ee , that is proportional to 
the neutrinoless double beta decay amplitude 46-48], 



i=l 



rmUt 



*2 



(50) 
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TABLE II. The three-flavour oscillation parameters within la 
error from Ref. [451 ]. 



parameters 



la 



lAmjjxl 
sin 2 6*i2 
sin 2 023 
sin 2 613 



(7.59to.S) x 10- 5 eV 2 



(2.40tg;i 2 ) x 10" 3 eV 2 



U. J18_ 016 



+0.07 



0.50 



< 0.035 at 90% C.L. 



the constraint from Tritium (3 decay [4Cj , m Vi , 

3 

m 2 Ve = rni\Uu\ 2 < (2.3eV) 2 at 95% C.L. , (51) 

2 = 1 

and the bound on the sum of light neutrino masses, T 
from cosmological and astrophysical data: 



If one refers to a more restrictive bound on T < 0.28 eV 
at 95% C.L. (53|, one then gets for the lightest neutrino 
mass the following restrictions: 



mi < 0.089 eV (NH) , m 3 < 0.084 eV (IH) . (55) 



The real and orthogonal matrix O e can be expressed 
in terms of the parameters e e a b defined in Eq. (1331) and 
the charged lepton masses are 



0.486661305 ± 0.000000056 MeV , (56a) 
102.728989 ± 0.000013 MeV , (56b) 
1746.28 ±0.16 MeV, (56c) 



T=J2 m *< °- 68 eV at 95% CL - • 



(52) 



This upper limit on T results from the combination of the 
Cosmic Microwave Background data of the WMAP ex- 
periment with supernovae data and data on galaxy clus- 
tering p9j . 

Bounds on \m ee \ can be estimated by taking into ac- 
count the best fit values of neutrino oscillation parame- 
ters and the upper limit of sin 2 #13 from Table |TT] together 
with the assumption of a hierarchical spectrum for the 
light neutrino masses [50 52], 



\m ee \ < 0.005 eV(NH) 



10 -2 eV < 



< 0.05eV(IH) . 



(53a) 



(53b) 



Since neutrino oscillations are only sensitive to the 
mass squared differences, the constraints given by 
Eqs. (I5T1) and (|52"|) can determine the lightest neutrino 
mass: mi in the case of normal hierarchy (NH) as in the 
charged lepton masses or m 3 in the case of inverted hier- 
archy (IH). The fact that these two hierarchies are still 
phenomenologically viable is due to the indetermination 
on the sign of the mass squared difference Am 2 ^. The 
constraint given by Eq. (1511) implies that mi (NH) or 771,3 
(IH) should be less than 2.3 eV, which is a rather poor 
constraint and needs to get an improvement in future 
experiments. 

The cosmological bound is much more severe and even 
if one takes the upper bound of quantity T given in 
Eq. (|52")l . one gets 



evaluated at Mz scale through the rcnormalisation group 
equations for QED in the MS scheme at 1-loop level |54l . 
|55| . In the limit where the parameters e e a b are small, the 
orthogonal matrix O e that diagonalises M e takes approx- 
imately [2^, [5^| the following form: 



(57a) 

(57b) 
(57c) 
(57d) 
(57e) 
(57f) 




mi < 0.22 eV (NH) , m 3 < 0.22 eV (IH) . (54) 



The unitary matrix, U v , that diagonalises the neutrino 
mass matrix mj given in Eq. (|46l) can be written in terms 
of the neutrino masses, the real parameter D v and the 
phase <p from Eq. (|41b|) . If one takes the phase ip = 
the matrix XJ V becomes an exact form of the remaining 
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parameters [lj|, with the matrix element moduli, 



{U v ) n 



{U v ) 



:>,:>, 



m,2mz{D v — m,\) 
D v (m 2 - mi)(m 3 - mi) ' 

mim3(m2 — D v ) 
D u (m 2 — mi)(m 3 — m 2 ) 

mim 2 [D v — m 3 ) 
£>y (ffi3 - mi)(m 3 - m 2 ) ' 

m 1 (m 1 - D u ) 
(to 2 — mi)(m 3 — mi) ' 

(D v -m 2 )m 2 
(m 2 - mi)(m 3 - m 2 ) ' 

m 3 (m 3 -D u ) 
(m 3 - mi)(m 3 - m 2 ) ' 

' mi(-D„ - m 2 )(D v - m 3 ) 
-Di/ (m 2 — mi)(m 3 — mi) ' 

m 2 (D v - mi)(m 3 - £>„) 
-Dy (m 2 - mi)(m 3 - m 2 ) ' 

m 3 (D v - mi)(D y - m 2 ) 
D v (m 3 - mi)(m 3 - to 2 ) ' 



(58a) 
(58b) 
(58c) 
(58d) 
(58e) 
(58f) 
(58g) 
(58h) 
(58i) 



and the signs for the matrix elements (U v )ij are in the 
case of normal hierarchy given by 



+ 



sign(mi) sign(m 2 ) sign(m 3 ) 
ysign(m 2 m 3 ) sign(mi) + 

and in the case of inverted hierarchy by 

/ + + + 

sign(mi) sign(m 2 ) sign(m 3 ) 
ysign(m 3 ) + sign(TOim 2 ) y 



(59) 



(60) 



In our numerics we have performed the full diagonalisa- 
tion with all possible values for tp by just applying the 
Eq. (HU). 

To study the new constraints that arise from a charged 
lcpton mass matrix M e in the NNI form and the ef- 
fective neutrino mass matrix m„ belonging to Class-II, 
we have varied all experimental charged lepton masses 
and neutrino mass differences within their allowed range 
given in Eq. (|5"6"|) and Table HH respectively. The mass 
of the lightest neutrino (mi in NH or m 3 in IH) was 
scanned for different magnitudes below 2 eV. In order to 
fully reconstruct the PMNS matrix, we have also varied 
the free parameters e e ab , D v and the phases Ki, k 2 , ip, 
defined in Eq. (|41[) . The remaining six parameters, 
A e , B e ,A v , B u , C e<v are determined from the values of 
lcpton masses and e e a b , D u , tp. The restriction in this 



scan was to accept only the input values which corre- 
spond to a reconstructed PMNS matrix U that naturally 
leads to the mixing angles £q 2 , # 23 and 6q 3 within their 
experimental bounds presented in Tabic |ITJ 





<<Sy4>: r 



0.002 0.004 0.006 0.008 
m.(eV) 



FIG. 2. Plot of sin 2 #13 as a function of mi in the case of 
Texture-II and normal hierarchy. The dashed lines correspond 
to sin 2 0i3 = 0.016 ± 0.010 at la C.L. from the global analy- 
sis [13] • The plot of sin 2 823 vs mi is not present here, since 
it reveals no correlation. 

From our search, we have found that the neutrino mass 
hierarchy can distinguish the two zero textures, II and 
H(i2), since Texture-II is only compatible with normal 
hierarchy, while TextureTI( 12 ) requires inverted hierar- 
chy. In Fig. [2j we plot sin 2 9i 3 against the lightest neu- 
trino mass m\ for the Texture-II where normal hierar- 
chy applies. For the mixing angles 6>i 2 and 9 23 we found 
no correlation at all. While in Fig. [3J we plot sin 623 , 
sin 2 #i 3 over the lightest neutrino mass, which is m 3 for 
the TextureTI( 12 ) since it is only compatible with in- 
verted hierarchy. In the plots of sin 2 9 13 over the lightest 
neutrino mass shown in Fig. [2] and [3] we have also drawn 
the new hint of non-zero sin 2 1 



'13 



'13 = 



0.016 ±0.010. 



(61) 



at la C.L. from the global analysis [57| of all avail- 
able neutrino oscillation data. A direct determination of 
sin 2 9i3 can affect significantly the validity of our model. 
It is indeed clear in Fig. 2] that neither Texture-II nor 
Texture-II( 12 ) are compatible with inverted or normal hi- 
erarchy, respectively. In both situations, the value of 
\Ui 3 \ is two orders of magnitude outside the required 
bound for sin 2 6*i 3 given in Table [LTl 

One can clearly see from Fig.[5]and Fig.[3]that for both 
textures the lightest neutrino mass, mi (NH) or m 3 (IH), 
is bounded. In the case of Texture-II one has, 



0.0013 eV < mi < 0.016 eV, 
whereas for Texture-II(i 2 ) one has, 

0.0042 eV < m 3 < 0.011 eV, 



(62) 
(63) 
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0.3 0.31 0.32 0.33 0.34 
sin 2 6„ 




J i I i I i L_ 



0.0O4 0.006 0.008 0.01 0.012 
m, (eV) 



Texture-II (NH) 



Texture-II (IH) 




Im l(eV) 




0.004 0.008 0.012 0.016 
m.(eV) 



0.004 0.006 0.008 0.01 0.012 
m(eV) 



FIG. 3. Plots of sin 623 versus sin #12 (left) and sin 813 
versus 7713 (right) in the case of Texture-II( 12 ) and inverted 



(12) 

hierarchy. The dashed lines correspond to sin 2 #13 
0.010 at lcr C.L. from the global analysis [57]]. 



0.016 ± 



Texture-II (IH) 



Texture-II,,, (NH) 



lite, 

"MM' 1 ' 




0.02 0.04 0.06 0.08 0.1 
m,(eV) 



0.02 0.04 0.06 0.08 0.1 
m,(eV) 



FIG. 4. Plots of sin 2 #13 as a function of the lightest neu- 
trino mass for Texture-II (left), 7713 for inverted hierarchy, and 
Texture-II(i2) (right), mi for normal hierarchy. The dashed- 
pointed line corresponds to the limit sin 2 #13 < 0.035 at 90% 
C.L. from the global analysis [451 ]. 



which let us to conclude that the light neutrino mass 
spectrum cannot be quasi-degenerated and there is no 
room to account for a massless neutrino state (as it is 
still allowed by general neutrino oscillation analyses). 
The upper bound on the constraint from Tritium /3 de- 
cay given in Eq. (|51[) does not have any impact in these 
results. Although the cosmological and astrophysical up- 
per bound on the sum of light neutrino masses seems 
rather severe, the upper limit on the lightest neutrino 
mass given in Eq. (|54[) or even the most restrictive value 
in Eq. (|55p are indeed above the bounds reported in 
Eqs. (jr ^l and (P I) . 



FIG. 5. Plotting |m ee | versus mi (left) and 771,3 (right) for 
Texture-II (NH) and -II(i2) (IH), respectively. 



Concerning the effective Majorana mass |m ee |, we 
present in Fig. [5] the value of |™ ee | as function of the 
lightest neutrino mass for both Texture-II and -H(i2) - 
From our scan, we have obtained the following limits for 

\m ee \ : 

6.4 x 1CT 4 eV < \m ee \ < 2.2 x 1CT 3 eV , (64) 

in the case of Texture-II and, 

0.015 eV < \m ee \ < 0.022 eV , (65) 

for Texture-II( 12 ) . These bounds obtained for |TO ee | are 
in full agreement with those given in Eq. (1531) . One ex- 
pect that future improvements on the experimental value 
of |m ee | may have impact on this textures (II and II(i2)) 
or on the Z4 symmetry itself, since a change on the mag- 
nitude of the experimental value of |TO ee | may drastically 
constrains our model. 

Before closing the section one may address the question 
which could be the smallest deviations to the Fritzsch 
ansatz, e e , for the charged lepton mass matrix M e ac- 
ceptable by the experimental data as it was done for 
the quark sector in Ref. (25|. In the quark sector, the 
lower bound for e q consistent with electroweak data can 
be evaluated by taking into account the values of the 
quark masses, CKM elements (V us , V c b , V u b) and the 
angle /3 = arg(— V c dV* b V t * d V t b) of the unitarity triangle 
listed in Table IIIII We have performed an update of the 
lower bound of the parameter e q defined in Eq. (|35[) , 



e„ > 0.188, 



(66) 



which has increased more than 10% compared with the 
value reported in [25[. 

The deviation from Hermiticity of the charged lepton 
mass matrix e e was computed for the Texture-II and - 



II 



(12) 



and the lower bounds e e > 0.0011 and e„ > 0.0013 
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TABLE III. Values of the quark masses, the CKM element 
moduli \V US \ , \V c b \ , \V u b\ an d the angle /3 of the unitarity tri- 
angle at the scale Mz ■ The quark masses are calculated at Mz 
scale by taking properly into account the 4-loop renormalisa- 
tion group equations for QCD in the MS scheme [H3.l55l.l5S' 
IriOl ] and using as input the masses given in [40j ] . 



1.4±0.5MeV 



0.62to.o? GeV 



m t = 170.2 ±1.0 GeV 



\Vus\ = 0.2253 ± 0.0007 



(41.0I^i) x 10' 



m d = 2.9 ±0.5 Me V 
m s = 58t\ 6 2 MeV 
m b = 2.86±g;S| GeV 

\V uh \ = (3.47±^ 6 2 ) x 10- 
sin 2/3 = 0.673 ±0.023 




TABLE IV. Summary of Texture-II and -II(i2) results. 

Texture-II (12) (IH) 

f -k * * 

★ 

V* o * 

sin 2 9is > 0.010 
0.0013 eV < mi < 0.016 eV 0.0042 eV < m 3 < 0.011 eV 
0.00064 eV < |m ee | < 0.0022 eV 0.015 eV < |m ee | < 0.022 eV 
e e > 0.0011 e e > 0.0013 



were found, respectively. These bounds imply that the 
neutrino data allow for the charged lepton mass matrices 
M e to be much closer to Hermiticity than the quark mass 
matrices, by two orders of magnitude higher. 



VI. CONCLUSIONS 

In this work we extended the flavour symmetry pro- 
posed in Ref. [2|| in the context of SU(5) Grand Unifica- 
tion. Such symmetry was imposed in the Lagrangian in 
order to force the quark mass matrices M u , Md to have 
NNI form, after spontaneous gauge symmetry breaking. 
In this GUT model, beyond the standard field content, 
three right-handed neutrinos and two Higgs quintets were 
added. It turns out that the light neutrino masses are 
generated through type-I seesaw mechanism and the low 
energy theory below the GUT scale is just a two Higgs 
doublet model. In this context, the minimal realisation 
of such discrete symmetry is Z 4 . 

Since in the SU(5) x Z4 model the charged lepton mass 
matrix, M e , are related to the down-quark mass matrix 
Md, it was not surprising to verify that the mass matrix 
M e acquired a NNI form, too. Instead, the right-handed 
neutrinos being singlets under SU(5), their Z 4 charges 
are free parameters and the neutrino mass matrix, m„ 
has no NNI form. There are in fact only six zero tex- 
tures allowed for the effective neutrino mass matrix, m v . 
However, the neutrino oscillation data select just two tex- 
tures among the six possibilities: Texture-II, which is 
compatible only with normal hierarchy for the neutrino 
mass spectrum and the Texture-II( 12 ), which demands 
neutrino mass spectrum to have inverted hierarchy. 

In Table IIV1 we summarise the main predictions for 
Texture-II and -H( 12 ). In both textures, the lightest neu- 
trino mass is predicted to be bounded and the neutrino 
mass spectrum to be hierarchical, therefore not compat- 
ible with a massless neutrino. Our results are also in 
agreement with three other relevant constraints, namely 
the effective Majorana mass m ee , the constraint from Tri- 
tium j3 decay and the cosmological bound on the sum 



of light neutrino masses. We have also obtained that 
the deviation of Hermiticity on the charged lepton mass 
matrices is two orders of magnitude lower than the one 
estimated for the quark sector. 

Future improvements on the knowledge of neutrino os- 
cillations, neutrinoless double beta decay, tritium beta 
decay and cosmological astrophysics measurements may 
be decisive for testing the viability of the SU(5) x Z4 
model. 
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Appendix A: Two-loop evolution of the running 
gauge couplings 

In this appendix we collect the two-loop renormalisa- 
tion group equations for the gauge coupling constants 



13 



oti (i — 1,2,3), which can be written in the form [6lM6 



The two- loop coefficients Cif are given by, 



d -1 h 1 , 

dt 1 2tt 8tt 2 ^ 3 3 

j 



(Al) 



fe=l,2 



where ai = 5/3 a y , 6^ are the usual one- loop beta co- 
efficients, bij and Cif are the two- loop beta coefficients. 
The quantities Tj 2 denote the quark and lepton Yukawa 
coupling matrices corresponding to the Higgs doublets 
$i,$ 2 - At the unification scale A, the gauge couplings 
Ui obey to the relation 



ai(A) = a 2 (A) = a 3 (A) . 



(A2) 



In the region of energy scales where one has only the 
SM degrees of freedom, the /3-function coefficients are 
given by: 
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C, 



if 
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3 3 1 

2 2 2 

\2 2 oj 



(A4) 



and they are neglected in our Runge-Kutta integration. 

Concerning the /3-function coefficients for the relevant 
particle content of the SU(5) theory, absent in the SM, 
one has the following coefficients for the doublets $1, $2, 
the colour triplets T\,Ti, the triplet S3 and the octet £§: 
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which are introduced at the appropriate intermediate 
scales. 
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